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ABSTRACT 



^ We have been able to compare with astrometric precision AstroDrizzle pro- 

cessed images of NGC 6720 (the Ring Nebula) made using two cameras on the 
Hubble Space Telescope. The time difference of the observations was 12.925 yrs. 
This large time-base allowed determination of tangential velocities of features 
within this classic planetary nebula. Individual features were measured in [N II] 
images as were the dark knots seen in silhouette against background nebular 
[O III] emission. An image magnification and matching technique was also used 
to test the accuracy of the usual assumption of homologous expansion. We found 
that homologous expansion does apply, but the rate of expansion is greater along 
the major axis of the nebula, which is intrinsically larger than the minor axis. 



We find that the dark knots expand more slowly that the nebular gas, that the 
distance to the nebula is 720 pc ±30%, and the dynamic age of the Ring Nebula is 
about 4000 yrs. The dynamic age is in agreement with the position of the central 
star on theoretical curves for stars collapsing from the peak of the Asymptotic 
Giant Branch to being white dwarfs. 

Subject headings: Planetary Nebulae: individual(Ring Nebula, NGC 6720)- 
instrumentation:miscellaneous:individual(HST,WFC3) 



1. Background and Introduction 

Using the spatial motion of the material in planetary nebulae (PNe) to understand 
their 3-D structure and evolutionary processes has been a goal of many studies of these 
ubiquitous objects. Olin C. Wilson (1950) first established the utility of radial velocity maps 
of a PN for determining its basic structure, showing that most objects could be explained 
as expanding ellipsoidal shells. Numerous studies since then have exploited this approach as 
high-velocity resolution spectrographs became more common and as high quantum efficiency 
CCD detectors replaced photograph emulsions. We now understand how the apparent form 
of many planetary nebulae (PNe) is the result interacting stellar winds. Slowly moving gas, 
expelled as the central star reaches the top of the Asymptotic Giant Branch, is overtaken 
by a fast wind of gas, expelled as the star begins to heat up. This is when the central star 
begins its collapse towards becoming a white dwarf. In the earliest phase the ejected material 
is mostly neutral and has a small ionized core. The shell may become fully ionized as the 
central star increases in temperature and remains at high bolometric luminosity. As the 
luminosity of the central star in ionizing photons drops, the outer parts of the material may 
again become neutral. From that point on the balance of ionized and neutral gas becomes a 
subtle balance. The decreasing gas density that results from the nebula's expansion causes 
the ionized fraction to increase. However, the decreasing luminosity of the star causes the 
ionized fraction to decrease. These factors mean that the dynamic structure that we see today 
is the result of multiple processes, rather than the simple ejection of a shell from a collapsing 
star. Nevertheless, it is worthwhile studying the dynamics of the PNe because it allows us 
to infer what must have been their earlier conditions and to predict their futures. 



^ased on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
NASA Contract No. NAS 5-26555. 



1.1. 3-D Structures from Radial Velocities 

Constraining determination of the 3-D structure and dynamics only to the use of line- 
of-sight (radial) velocities is an important limitation. Interpretation of the observed radial 
velocities usually demands a simplifying assumption about the general nature of the motion. 
Most often the assumption is of homologous expansion; that is, that the spatial velocity away 
from the central star scales with the spatial distance. Homologous expansion is like explosive 
expulsion of particles of varying velocity. Knowing now that the expanding material is shaped 
and accelerated by multiple processes makes it surprising that PNe resemble homologous 
expansion, a result going back to Wilson's (1950) early study. Wilson showed that emission 
originating farther from the central star shows a higher velocity of expansion. Mapping 
the radial velocities of PNe has become straightforward and common, but these data only 
map one component of the spatial motion (Vr). However, one cannot underestimate the 
importance of Wilson's (1950) study since quantitative efforts to determine the 3-D structure 
of PNe originate with this paper. 



1.2. Determining Tangential Velocities in Planetary Nebulae 

Potentially one can create better models for the PNe if we also have measurements of 
their tangential motions. However, for a fixed value of the spatial motion perpendicular 
to the line-of-sight (Vr) the apparent motion in the plane of the sky (n) decreases with 
increasing distance (D) to the target, and many fewer nearby candidate objects exist. When 
expressed in the common units of Vr in kilometers/sec, // in "/yr, and D in parsecs, the 
relation is /i = 0.211 xV^/ D. Inserting representative values of Vr = 20 km s _1 and D = 
500 pc gives an expected /i of 0.0089 "/ yr , 8.9 milliarcsec (mas) per year. Determination of 
H was first attempted by Liller et al. (1966) using ground-based images of planetary nebulae 
obtained with several large telescopes over a period of up to 62 years. No believable values 
were produced because of the blurring by astronomical seeing (characteristically about 2" 
in the older images employed) and the lack of a clear filtering of single emission-lines. The 
latter limitation is because different emission-lines can arise from different positions in the 
nebula, thus blurring the image. 

As the Very Large Array (VLA) radio interferometer began to produce stable images 
of 1" resolution and better, the possibility was recognized that this telescope could provide 



usable /z values by special processing of the data (Masson 1986 Gomez et al. 1993 Hajian 



et al. 1993, 1995). A summary of results was given by Terzian (1997). Limited to the study 
of high surface-brightness PNe and unable to give unambiguous results for different portions 
of the PNe, the method is no longer in wide use. 



The VLA method of obtaining \x has essentially been succeeded by the use of two-epoch 
images using the Hubble Space Telescope (HST). With its stable optical- wavelength images 
of slightly better than 0.1" resolution and its ability to clearly isolate regions emitting specific 
emission-lines, this becomes the more promising approach. Reed et al. (1999) employed pairs 
of images of NGC 6543 using the high resolution CCD camera (45.5 mas pixels) of the Wide 
Field Planetary Camera 2 (WFPC2) with a time separation of 2.92 yrs. They identified 
three approaches. The first approach was to measure radial profiles of samples on the two 
sets of images and then look for displacements using a gradient recognition algorithm. The 
second approach also used radial profiles, but fit gaussian profiles to features and looked for 
the changes in their centers. The final technique was to magnify the first epoch images by a 
succession of scale factors, align their centers on the central star, and recognize at what scale 
factor the magnified first epoch and the second epoch images matched. This last approach 
utilizes all of the nebula's image, but its accuracy is difficult to quantify and works only 
in the case of near-homologous expansion. Motions of up to 6.5 mas yr _1 were found. In 



a subsequent study (Palen et al. 2002) of four PNe with time differences of 3.9 yrs to 4.2 
yrs, the first and second approach were used to successfully determine \i in three objects. 
However, they compared broad bandpass images (the WFPC2 F555W) with monochromatic 
images (the WFPC2 F502N), arguing that the signal of both filters was dominated by the 
[O III] emission line at 500.7 nm. Values up to 2 mas yr _1 were found. 



1.3. Tangential Velocities in the Ring Nebula 

As the archetypical elliptical-form PN, NGC 6720 (henceforth the Ring Nebula) has 
attracted considerable attention in observational programs. However, it is significantly lower 
in surface-brightness than most of the PNe that have been the subject of radial velocity and 
tangential motion studies, and high resolution spectroscopy has been a relatively recent 
subject. O'Dell et al. (2002) used WFPC2 [N II] emission-line images with a 1.944 yr time 
difference and the magnification technique of Reed et al. (1999) to determine a magnification 
factor (M) of M=1.00067±0.0001 yr -1 . This study employed the coarser (100 mas) pixels of 
the wider field of view CCD's of the WFPC2. With an average diameter of its Main Ring 
of 36", this corresponds to an expansion of about 24 mas yr -1 . 

In a discussion of possibly homologous expansion it is more convenient to express the 
results in terms of what we call here the expansion scale factor E, where E = fi/<p, where <fi is 
the angular distance from a point in the nebula to the central star. If homologous expansion 
applies, then E will be constant. It is analogous to the Hubble constant that used to describe 
to the first order of accuracy the expansion of the universe. The expansion scale factor is 



most conveniently expressed with units of mas yr 1 / // . For the magnification technique 
E=(M-l)xl0 3 . The O'Dell et al. (2002) results can then be summarized as E=0.67±0.1 
mas yr" 



-y. 



However, a study of radial velocities in the Ring Nebula (O'Dell et al. 2007) questioned 



the accuracy of this result because the second epoch images had the nebular centered on 
the central star, which meant that the image of the Main Ring fell across all four CCDs in 
the WFPC2, whereas the first epoch images fell onto only a single CCD. Since the relative 
positions of the CCDs are more difficult to determine than the astrometric corrections within 
a single CCD, there could be a systematic error. This limitation was circumvented in a later 



study (O'Dell et al. 2009), where single CCD images with a time difference of 9.557 yr 



were employed. All of the techniques introduced by Reed et al. (1999) were used with the 
result that E=0.23±0.1 mas yr -1 /" for [N II] emission and four dark features seen in F502N 
gave E=0.30±0.17 mas yr -1 /". These results indicate that the O'Dell et al. (2002) numbers 
were inaccurate, probably because of using multiple CCDs in the second epoch images. The 
O'Dell et al. (2009) result also argues for the first time that the dark knots that mark the 
Ring Nebula and are ubiquitous in nearby PNe ( O'Dell et al.|2002 ) expand at approximately 
the same rate as the nebular material. The big uncertainty in the values of E means that a 
conclusion about commonality of motion is quite uncertain. 



1.4. 3-D Modeling and the Value of Tangential Velocities in the Ring Nebula 

Attempts to model the Ring Nebula in 3-D have a long history. These are summarized 



in a recent publication (O'Dell et al. 2013a) that draws on the extensive imaging and radial 



velocity mapping that has been done. In that study, it was found that the Main Ring of the 
nebula, where most of its optical light originates, is a non-circular disk of ionization bounded 
gas. Near the central axis of the disk material extends towards the observer in the form of 
lobes, seen only in projection onto the plane of the sky. The Main Ring is surrounded by 
a glow of [O III] emission and two low ionization halos. All three of these last features 
are visible as fossil radiation, material that was previously photoionized but is now shielded 
from the ionizing radiation of the central star. Hundreds of dark knots are best seen in 
extinction against the Main Ring's [O III] emission with the F502N filter, although many of 
them have low ionization arcs on the side of the knot that faces the central star. The dark 
knots generally point toward the central star, but do not have the radial symmetry of form 



seen in the Helix Nebula (O'Dell et al. 2005). 



Since the polar axis of the Ring Nebula is pointed within about 10° of the observer, 
the Main Ring is viewed nearly in the plane of the sky and tangential velocities can be 



determined. These tangential velocities can be used to understand the structure and distance 
of the nebula. The Ring Nebula is also a prime target for addressing the important question 
of the relative motion of the nebular gas and the dark knots. 

As part of HST program GO 12309 high signal-to-noise ratio images were made in 
the strongest emission lines of the entire nebula using the Wide Field Camera 3 (WFC3) 
that replaced the WFPC2 during the last servicing mission. These images were used to 
create the 3-D model discussed in O'Dell et al. (2013a). Images were also with smaller 
field-of-view diagnostic filters, and these were used to determine the physical conditions 
(ionization, temperatures, densities) within the nebula ( |Q'Dell et al. 2013b). The WFC3 



pipeline-processed images sample better (40 mas per pixel) the point-spread-function of 
the HST and thus provide the highest resolution images of the Ring Nebula. As a new 
camera, no earlier observations exist of the Ring Nebula. However, one can accurately 
use the earlier WFPC2 images as first epoch images because the pipeline-processed images 
extracted from the HST archives have had an accurate distortion-correction applied, thus 
rendering the early and recent images geometrically "flat". The geometric corrections step 
is called "AstroDrizzle" and is applied to both the old and new images. It is only necessary 
to scale the WFPC2 images to the same pixel size as the WFC3 images, align the central 
star images, and rotate them into alignment. One then has astrometrically useful paired 
images whose alignment accuracy should be about 0.2 WFC3 pixels. Not only are these 
recent images of higher spatial resolution, they offer the longer time base of 12.925 years 
when combined with the first images made in program GO 7632. 



In this paper we describe the method of aligning and measuring the image pairs (§ 2.1), 
present the results for measurements of gaseous features in [N II] and the dark knots in 
F502N (§ 2.2 and 2.3), and discuss what this tells us about the Ring Nebula (§[3]). 



2. The Images 

The images used in this study were obtained in the HST programs GO 7632 (WFPC2, 
1998-10-16) and GO 12309 (WFC3, 2011-09-19). Multiple exposures were made at each 
epoch, thus allowing cosmic ray event correction. The images in filter F502N isolate well 
the 500.7 nm line of [O III], and the images in filter F658N isolate well the 658.3 nm line of 
[N II] . The [N II] emission arises from a thin layer close to the main ionization front of the 
nebula (O'Dell et al.|2013a). In contrast, the [O III] emission arises from an extended region 



lying between the [N II] emitting layer and the Hell core in the center of the nebula (O'Dell 



et al. 2013a). This means that the F658N images provide a means of tracing the motion of 



gas near the main ionization front, while the primary value of the F502N images lies in their 



ability to show the dark knots in detail. 



2.1. Combining the WPC2 and WFC3 Images 

The first epoch GO 7632 AstroDrizzle images were downloaded from the HST data 
archive and put into the same pixel scale as the second epoch GO 12309 AstroDrizzled by 
using the IRAF q task "magnify" and the scaling factor of 2.5262758. The scaling factor 
could not be determined from the few stars in both Ring Nebula images was obtained by 
comparing rich star fields in the Orion Nebula Cluster imaged with the WFPC2 in program 
GO 5085 and with the WFC3 in program GO 12543. The two sets of images were aligned 
with IRAF task "geomap/rotate" using the four stars available. The alignment accuracy 
of the central star is about 3 mas, and the determination accuracy of the rotation is about 
0.015°. 

The aligned images were scaled to the same signal level, and the WFC3 image was 
blurred to the same stellar image point-spread-function as the WFPC2 image using the 
IRAF task "gauss" . A ratio of the first and second epoch images was made, with the second 
epoch as the denominator. These ratio images revealed motions that had occurred during 
the 12.925 yr interval between the observations. For bright objects, the F658N ratio image 
shows motion of objects in the direction of the dark edge. For the dark knots, the ratio 
image in F502N shows motion of objects in the direction of the bright edge. These images 
are shown in Figure [TJ 



2.2. Measuring the Motion of Features with a Least-Squares Method 

The ratio images were used as a guide in determining samples to be measured using 
a least-squares code developed by Hartigan et al. (2001) from an approach originated by 
Currie et al. (1996) that compares incrementally shifted images. The motion is assumed to 
be correct when there is a minimum difference between the reference and shifted image. In 
this paper we call this method the ZSQ method. The results for both small dark features 
seen in the F502N images and large bright features seen in the F658N images are shown 
in Figure [2} Results for individual determinations of /i are given in Table 1 and Table 2. 



2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associ- 
ation of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science 
foundation. 



After recognizing patterns in the tangential velocities we grouped the \i values into quadrants 
within the Main Ring as shown in Figure [2] The averages within the quadrants are presented 
in Table 3. 

It is possible to approximate the uncertainties of the measurements. The ZSQ gives 
much smaller uncertainties than indicated by a comparison of \x for two nearby objects. 
Excluding the F502N dark features in the SW quadrant (c.f. §[3|, the average motion of all 
features is 2.12±0.33 ACS pixels. Over the time interval of the observations, this corresponds 
to ft = 6.56 mas yr _1 . If the alignment uncertainty is that of the central star image ( 3 mas 
or 0.08 pixels), that component of the uncertainty of the magnitude of \x is only 0.2 mas yr _1 
and an uncertainty in the direction of motion of 2.2°. The average distance of this set of 
samples is 27.5" (688 pixels). The uncertainty in the rotation between the two sets of images 
(0.015°) would cause the tips of the features to be displaced by 0.18 pixels, which would 
produce an error in the direction of motion of 4.9°. These two considerations mean that 
uncertainty in the alignment propagates into small uncertainties in the derived magnitude 
and direction of the motion of features in the nebula. Typically, the uncertainties in the ZSQ 
solutions are about 0.05 pixels, which corresponds to a 3% uncertainty in the magnitude of // 
and an uncertainty of about 1.5°. The significant differences in the magnitude and direction 
of \i for nearby objects in the Main Ring are likely to be real and can represent significant 
small-scale fluctuations in the gas flow of the Main Ring. An alternative explanation may 
be that two adjacent lines of sight are dominated by radiation from very different parts of 
the 3-D Main Ring. 



2.3. Measuring Motions with the Magnification Method 

We adopted a variation of the magnification method of determining motions. The un- 
derlying assumption of the method is that the tangential motions are homologous. When 
a correctly magnified first epoch image is aligned with the second epoch image (using the 
central star) and a ratio of the two is made, a flat image of constant value should result. 
However, if the expansion scale factor is not constant in all parts of the nebula, the em- 
ploying a single optimum magnification (where the ratio of images becomes constant) will 



not represent fully what is happening. Recall from § 1.3 that E and the magnification (M) 
are related by E=(M-l)xl0 3 when E is conveniently expressed in units of mas yr -1 /". We 
determined the optimum magnification and E for individual regions within the image of the 
Main Ring, as shown in Figure [3j The F502N samples were selected from where the nebular 
background varied little and contained numerous dark knots. The results should represent 
the motion of the dark knots. The F658N samples were selected from regions showing struc- 



ture in the [N II] image. The results should indicate motion of ionization fronts aligned with 
the observer's line-of- sight. 

Several steps were used to derive the optimum magnification. Initially, the first and 
second epoch images were scaled to the same signal level. The first epoch images were 
magnified in steps of 0.00005. The central star was aligned on the resultant images, and 
a ratio of the two images was made for each magnification. Finally, the IRAF task "im- 
stat" was applied to establish the standard deviation (STDDEV) of the signal within the 
sampled region. There was always a difference between the two images. The STDDEV was 
largest when the magnification was furthest from optimum, because the ratio-image was the 
result of comparing two distinct images, each with its own characteristics. However, as the 
magnification factor reached the optimum value, the STDDEV reached a minimum, then 
increased at values larger than the optimum magnification. The STDDEV was never zero 
because the signal-to-noise ratio of the first and second epoch images was not the same. 
Imperfect matching of the image resolution would also contribute to a non-zero STDDEV 
value. The images had been approximately matched to the same point-spread-function, but 
subtle differences in fine-structure would exist and contribute to a non-zero STDDEV even 
at the optimum magnification. The results of the optimum magnification and the expansion 
scale factor determinations are given in Table 4 and Table 5. The group averages of the 
magnification results are given in Table 6. 



3. Discussion 

3.1. Large Scale Properties of the Ring Nebula's Main Ring 

The results of the analysis of the motions was made in terms of average values of E. 
The quadrants we adopt are primarily meant to recognize regions of different motion, but 



also they agree with the natural divisions suggested by the most current 3-D model (O'Dell 



et al. 2013a). A remarkable symmetry of features appears in opposite quadrants, which 
indicates similarity of properties. This justifies grouping the results into the quadrants along 
the major (NE+SW) and minor (NW+SE) axes. We adopt as the average of the ZSQ and 
magnification methods <E([N II], minor) > = 0.22±0.04, <E([N II], major) > = 0.27±0.05, 
<E(F502N,minor)> = 0.20±0.03, and <E(F502N,major)> = 0.18±0.02. 

Our measurements of \x are limited to the Main Ring of the nebula and include no 
measurements of inner halo or outer halo features. O'Dell et al. (2007) first established and 
O'Dell et al. (2013a) confirmed the ellipticity of the Main Ring. The nebula's equatorial 
concentration of material is elliptical in form, rather than the apparent ellipticity being the 



10 



result of a circular ring seen at an angle. The small difference of the major and minor axes 
<E([N II] )> values is consistent with the intrinsically elliptical equatorial concentration of 
the Main Ring. The minor axis material is expanding at a slower rate either because it 
was subject to less force during the phase when the Main Ring was being formed, or the 
over-burden of outer neutral material is greater there. 



3.2. Relative Velocities of the Dark Knots and the Nebular Gas 



Systematic differences of the averaged E values in [N II] and F502N are seen to exist 
along both axes. This difference is greatest along the major axis, where <E([N II] )> = 
0.27±0.05 and <E(F502N)> = 0.18±0.02. The difference along the minor axis is much 
more uncertain, with <E([N II])> = 0.22±0.04 and <E(F502N)> = 0.20±0.03. A smaller 
value of E(F502N, dark knots) as compared with E([N II], bright gas) indicates that the 
dark knots do not expand as fast as the ionized gas in which they exist. 

Two families of models are used to explain the origin of the dark knots found in PNe. 



The first (Dyson et al. 1989) argues that the knots are condensations of dust and gas that 



originated in the extended atmosphere of the central star while it was still at the top of 
the Asymptotic Giant Branch. The similarity of the E values found today argue against 
this origin. It would be highly unlikely that knots sufficiently dense to survive the first 
phases of shell ejection would receive a radial acceleration this similar to that of the low 



density gas. The second set of models (Capriotti 1973 Vishniac 1983 Garcia-Segura et 



al. 2006 ) invoke formation of the knots much later in the life of the nebula. Considerable 



observational evidence (O'Dell et al. 2002 O'Dell et al. 2004) indicates that the knots arise 



near the main ionization front. Our results show that the knots, once formed, appear to lag 
behind the motion of the surrounding gas, and this feature must be factored into theories of 
their formation and calculations of the subsequent properties. 

Previously Meaburn et al. (1998) argued that the system of well-defined dark knots in 
the Helix Nebula are expanding at a rate of about one-half that of the nebular gas. That 
conclusion was drawn from study of the radial velocities of the small crescents of ionized gas 
formed on the side of the knots facing the central star. The primary concentration of gas in 



the Helix Nebula is in an equatorial ring lying almost in the plane of the sky (Meaburn et 



al. 1998; O'Dell et al. 2004), much like the Ring Nebula. Expansion motions are nearly in 



the plane of the sky, meaning that results determined from radial velocities are inherently 
less valuable. A study of the tangential velocities should be more accurate for addressing the 
question of the relative motion of the knots. However, it appears that the study of tangential 
velocities of the Ring Nebula and the radial velocities in the Helix Nebula come to the same 
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conclusion: the knots are expanding more slowly and are decoupled from the gas motion. 



3.3. The Distance to the Ring Nebula 



A distance (D) for the Ring Nebula can be obtained from an accurate knowledge of 
its dynamical properties. For homologous expansion, the relation that applies is D(pc) = 
(ki x 1Z x cos#)/(4.74 E), where ki is the scale factor relating the spatial expansion velocity 
(km s _1 ) to the true (not projected) distance to the central star (expressed as "), 1Z is the 
pattern velocity (the ratio of the ionization front velocity compared with the gas velocity), 9 
is the angle of the equatorial plane of the nebula to the sky (°), and E is the expansion scale 
ratio. Adopting the values in O'Dell et al. (2007) of ki = 0.65 (km s" 1 / ") and K = 1.34, 
the relation becomes D(pc) = 179xE _1 , where E is in mas yr -1 /". From this work we adopt 
the average of the major and minor axes for [N II] (E=0.25 mas yr -1 /") and find D = 716 
pc. It is difficult to identify an uncertainty in this distance because it is dependent upon the 
validity of the dynamical model used. The key uncertainties are in the determination of ki 



(O'Dell et al. 2007) and TZ (O'Dell et al. 2009). The resulting uncertainty in the distance 



is probably about 30%. The derived distance should be compared with the astrometric 
distance of 700±2qo pc ( Harris et al.|2007 ). Our rounded value of 720 pc is the best distance 
to be adopted. This is only a small change from the previously recommended value of 740 



pc (O'Dell et al. 2009). 



3.4. Timescales 



The reciprocal of the expansion ratio has the dimension of years and represents the time 
since the onset of the assumed homologous expansion. Adopting the average of E = 0.25 



£3.3) gives 4000 yrs since onset of the expansion. This method of calculation produces an 



overly simplified dynamical age and only has the meaning as a measure of when the basic 
dynamics properties of the expanding nebular shell were established under the assumption 
that the size at that point was small. 

The current characteristics of the central star are (log L/L ) = 2.3±0.2 and log T « 
5.1±0.02) ( O'Dell et al.|200"7 ) and comparison with the theoretical path for a collapsing star 
of about M = 0.62 M (B16cker 1995) indicates that the time since onset of stellar collapse 



is about 5000 yrs, i.e. similar to the dynamic age. 
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4. Summary and Conclusions 

We have been able to utilize early HST WFPC2 images of the Ring Nebula in combina- 
tion with recent WFC3 images to determine tangential motions of [N II] emission features 
near the main ionization front of this ionization-bounded nebula. We could also determine 
the tangential velocities of the dark knots seen as extinction of the diffuse nebular [O III] 
emission. 

Well defined radial expansion patterns were found in both the [N II] emission-line fea- 
tures and the dark knots. The results from a second method of determination of the motions 
under the assumption of homologous expansion are in excellent agreement with the results 
from measurement of individual features. 

The expansion scale ratio E that characterizes the scale of motions within the nebula 
is not the same on both axes of the Main Ring. It is larger for the major axis, a result 
consistent with the inherent ellipticity of the Main Ring. The apparent ellipticity of the 
Main Ring is primarily due to its inherent shape, rather than being the result of a circular 
object being observed at an angle to the plane of the sky. The suggestions of spiral structure 



are the results of the asymmetries of the boundaries of the Main Ring and the lobes (O'Dell 



et al. 2013a). 



The reciprocal of the expansion scale factor indicates that the dynamic age of the Ring 
Nebula is about 4000 yrs. Combining the new results for the expansion factor and the previ- 
ously established spatial expansion velocity scale factor determined from radial velocities has 
allowed the calculation of a distance to the Ring Nebula of 720 pc, with a likely uncertainty 
of 30 % . 
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Fig. 1. — These two images are both 120" by 104" and have been cropped to show only the 
region common to both the first epoch and second epoch images. They show the ratio of 
the first epoch divided by the second epoch images in the F502N and F658N filters. The 
method of their derivation is described in § [2] For bright objects, the F658N ratio image 
shows motion of objects in the direction of the dark edge. For the dark knots, the ratio image 
in F502N shows motion of objects in the direction of the bright edge. The nearly horizontal 
double band of small dark dots across the middle of the image are caused by cosmic ray 
artifacts in the region of the WFC3 image for which adequate cosmic ray artifact correction 
could not be made. The location of the central star is indicated by a black circle. Note 
that the high proper motion star enclosed in the white oval was not used in aligning the two 
epoch images. 
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Fig. 2. — The same fields of view as in Figure [T] are shown. The motions determined by the 
ZSQ method described in § [2] are shown as vectors. At an assumed distance of 740 pc, a 
motion of one ACS pixel corresponds to a tangential velocity of 10.9 km s" 1 . The velocity 
scale for the vectors is shown on both panels. The small samples used for determination 
of the dark knots' motions are shown as open circles in the F502N image dominated by 
background [O III] emission and the larger samples used for determination of the bright 
[N II] features' motions are shown as open boxes in the [N II] image. 




Fig. 3. — The same fields of view as in Figure [T] are shown. The areas used in determining 
the optimum magnification are designated by heavy closed rectangles. 
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Table 1. ZSQ Method Motions for [N II] Bright Features 

Sample-PA* (°) Sample-Distance (") \x (mas yr" 1 ) Vcctor-PA* (°) E (mas yr" 1 /") Quadrant 

021.9 30.8 6.51 037 0.211 NW 

007.1 28.9 

344.8 25.7 

339.7 28.8 
328.1 24.8 

327.1 27.8 
320.4 28.2 

287.2 27.4 

277.4 35.1 

246.5 29.3 ,.12 2(17 0.213 S\V 
236.5 40.7 
225.5 30.9 

208.0 29.4 

195.1 21.5 

194.3 31.5 5 71 1.9(1 0.I.M Sli 
187.3 31.8 
173.5 34.1 

161.8 26.7 

150.9 27.2 
127.3 32.0 
114.0 23.2 

099.0 29.5 10. D 1.09 0.353 NE 

083.5 39.9 

069.9 33.9 

067.3 26.1 

042.5 25.3 



Note. — *PA means the Position Angle, measured counter-clockwise from north. 



6.51 


037 


0.211 


5.11 


027 


0.177 


6.66 


352 


0.259 


7.97 


359 


0.277 


7.36 


342 


0.297 


4.14 


026 


0.167 


6.08 


340 


0.216 


5.85 


305 


0.213 


5.75 


297 


0.164 


7.12 


267 


0.243 


11.63 


239 


0.286 


7.63 


210 


0.247 


6.88 


209 


0.234 


7.25 


198 


0.338 


5.71 


196 


0.181 


5.79 


195 


0.182 


5.35 


186 


0.157 


7.01 


159 


0.263 


5.91 


161 


0.217 


6.03 


126 


0.189 


7.17 


111 


0.309 


10.41 


109 


0.353 


7.34 


127 


0.184 


7.29 


075 


0.215 


7.00 


059 


0.268 


7.33 


035 


0.290 



Table 2. ZSQ Method Motions for F502N Dark Features 



Samplc-PA* (°) Sample-Distance (") 



014.4 
007.8 
002.4 
334.3 
310.0 
300.9 
296.7 
280.1 
274.1 
266.3 
261.6 
252.6 
248.4 
240.0 
233.6 
233.1 
227.6 
219.2 
209.6 
205.5 
201.3 
198.7 
183.2 
168.7 
148.0 
153.1 
129.8 
104.0 
098.0 
073.5 
062.5 
048.6 
034.9 



16.3 
24.2 
25.5 
25.8 
25.4 
21.8 
28.0 
24.1 
21.1 
25.7 
22.4 
28.1 
25.1 
30.6 
28.0 
34.8 
32.8 
28.6 
27.5 
26.7 
27.9 
27.5 
25.7 
26.1 
25.0 
25.3 
28.6 
27.8 
29.8 
26.1 
23.3 
24.8 
21.2 



as yr ) 


Vcctor-PA* (° 


) E (mi 


5.45 


045 


0.334 


7.16 


010 


0.295 


5.86 


032 


0.230 


6.63 


013 


0.257 


4.92 


010 


0.194 


5.63 


002 


0.258 


5.23 


336 


0.188 


4.01 


329 


0.166 


3.40 


013 


0.161 


2.73 


331 


0.089 


2.54 


347 


0.113 


1.64 


247 


0.058 


2.81 


296 


0.112 


1.26 


323 


0.041 


1.80 


336 


0.064 


2.05 


271 


0.059 


2.61 


237 


0.080 


1.18 


202 


0.041 


3.28 


222 


0.119 


1.76 


169 


0.066 


0.93 


093 


0.033 


1.34 


173 


0.049 


6.79 


145 


0.264 


6.77 


163 


0.253 


4.75 


123 


0.190 


5.48 


126 


0.217 


4.91 


135 


0.172 


5.30 


098 


0.191 


4.77 


110 


0.160 


9.09 


095 


0.348 


7.26 


076 


0.312 


6.81 


066 


0.274 


6.96 


053 


0.329 



V") Quadrant 



NW 



SW 



SE 



NE 



Note. 



*PA means the Position Angle, measured counter-clockwise from north. 
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Table 3. Averaged ZSQ Motions Grouped by Quadrants 



Region 



Average fi (mas yr 1 ) Average Expansion Scale Factor (mas yr 1 /") 



[N II] Bright NW 


6.2+1.2 


[N II] Bright SW 


8.0+2.0 


[N II] Bright SE 


6.1+0.7 


[N II] Bright NE 


7.9+1.4 


[N II] Bright NW+SE 


6.2+0.9 


[N II] Bright NE+SW 


8.0+1.6 


F502N Dark NW 


4.4+1.2 


F502N Dark SW 


4.3+1.1 


F502N Dark SE 


5.7+1.6 


F502N Dark NE 


4.1+1.1 


F502N Dark NW+SE 


4.9+1.4 


F502N Dark NE+SW 


4.2+1.1 



0.22+0.05 
0.27+0.04 
0.21+0.05 
0.26+0.07 
0.22+0.05 
0.27+0.05 
0.19+0.04 
0.15+0.04 
0.22+0.06 
0.17+0.05 
0.20+0.05 
0.16+0.04 



Table 4. Magnification Method Results for [N II] Bright Samples 



Sample 


Optimum Magnifier 


E (mas yr 1 /") 


Quadrant 


a 


1.00275+0.00025 


0.21±0.02 


NW 


b 


1.0035 ±0.00025 


0.27±0.02 


NW 


c 


1.00288+0.00025 


0.22±0.02 


NW 


d 


1.00275+0.00025 


0.21±0.02 


SW 


c 


1.00425+0.00025 


0.33±0.02 


SW 


f 


1.0025 ±0.00025 


0.19±0.02 


SE 


g 


1.00325+0.00025 


0.25±0.02 


SE 


h 


1.0025 ±0.00025 


0.19±0.02 


SE 


i 


1.00425±0.00025 


0.33±0.02 


NE 


J 


1.00325±0.00025 


0.25±0.02 


NE 



Note. — *Sample regions are shown in FigurepS] 
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Table 5. Magnification Method Results for F502N Samples 



Sample 


Optimum Magnifier 


E (mas yr" 1 /") 


Quadrant 


A 


1.0025 ±0.00025 


0.19±0.02 


NW 


B 


1.00275 ±0.00025 


0.21±0.02 


NW 


C 


1.00225±0.00025 


0.17±0.02 


NW 


D 


1.00225±0.00025 


0.17±0.02 


sw 


E 


1.00225±0.00025 


0.17±0.02 


SW 


F 


1.0025 ±0.00025 


0.19±0.02 


SW 


G 


1.00275±0.00025 


0.21±0.02 


SE 


H 


1.00275 ±0.00025 


0.21±0.02 


NE 


I 


1.0025 ±0.00025 


0.19±0.02 


NE 


J 


1.00225±0.00025 


0.17±0.02 


NE 



Note. 



* Sample regions are shown in Figure [3] 



Table 6. Averaged Magnification Method Results Grouped by Quadrants 



Region 



1-Average Optimum Magnifier Average Expansion Scale Factor (mas yr /") 



[N II] Bright NW 


3.04±0.40 


[N II] Bright SW 


3.50±0.70 


[N II] Bright SE 


2.80±0.35 


[N II] Bright NE 


3.75±0.5 


[N II] Bright NW±SE 


2.73±0.30 


[N II] Bright NE+SW 


3.62±0.75 


F502N Dark NW 


2.50±0.25 


F502N Dark SW 


2.33±0.14 


F502N Dark SE 


2.75±0.25 


F502N Dark NE 


2.5±0.25 


F502N Dark NW+SE 


2.60±0.22 


F502N Dark NE+SW 


2.42±0.20 



0.23±0.03 
0.27±0.06 
0.21±0.04 
0.29±0.04 
0.22±0.03 
0.28±0.06 
0.19±0.02 
0.18±0.01 
0.21±0.02 
0.19±0.02 
0.20±0.02 
0.18±0.02 



